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Chapter I 


INTRODUCTION 


Since the advent of a systematic procedure for designing 
automatic control systems a generation ago, automatic control 
systems are becoming more common. At this time one can 
foresee an evolutionary step whose importance is comparable 
with that of a generation ago. This important step is the 
introduction of digital computers to replace analog computing 
elements in these systems. With the advancement of digital 
technology, digital computers have continuously become more 
efficient, expanding into new applications with each major 
technological improvement. Recently the power of digital 
computers has been made available at the integrated circuit 
level. They compute very fast, and are multifunctional, 
small and reliable. They give rise to the possibility of 
using dedicated digital processing units for a variety of 
complicated tasks; for example, the inclusion of digital 
processors as parts of various process control systems. 

This report presents the results of a study which in- 
dicates that replacing the analogue control system for the 
MIT-NASA prototype magnetic balance with a digital system 
is feasible. Preliminary tests with a one- dimensional 
suspension system were successful and it appears that a 
digital control system can greatly improve the speed, 
accuracy and versatility of magnetic balance operation. 

Magnetic balance and suspension systems for wind 
tunnels have been in use for over a decade (1,2) and have 
proved useful not only because of their accuracy but also 
because they offer the opportunity to conduct classes of 
tests that have been difficult to carry out in the past. 

The basic idea of the MMSB can be explained using the 
functional block diagram in Figure 1. As it is shown, the 
suspension system is a position control servo, designed to 
maintain the suspended model at some set position within 
the suspension magnet structure, according to the position 
command signals. The position of the model is translated 
into electrical signals by means of a position sensing de- 
vice, and these signals are compared with the position 
command signals. The resulting position error signals are 
modified in passing through a set of feedback compensation - 
networks . 
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The compensated error signals are then amplified and 
applied to the suspension system magnet coils. The magnetic 
fields produce magnetic forces on the suspended model accord- 
ing to the "magnetic force and moment relations" which are 
functions of the model position and orientation, as indicated 
by the internal feedback path. The magnetic forces and 
moments added to the aerodynamic and gravity forces and 
moment, act through the model inertias to cause changes in 
the position and orientation of the model. These changes are 
in turn sensed by the position sensor system and the servo 
loop is closed. The function of the servo loop is to minimize 
the integral of position error by continuously counteracting 
the aerodynamic and gravity loads. 

In addition to maintaining model position it is necessary 
to determine the aerodynamic loads acting on the suspended 
model. As we can see, the unknown aerodynamic forces and 
moments applied to the model by the wind are balanced against 
gravitational forces and moments, and magnetic forces and 
moments. Since gravitational forces and moments are known, 
the magnetic forces and moments can be calculated from elec- 
trical currents flowing through the magnetic coils and 
appropriate calibrations. As a result the aerodynamic forces 
and moments are determined within the limits of experimental 
accuracy , 

The following sections of the report describe: 

1. The present analog control MMSB system at the 
M.I.T. Aerophysics Laboratory. 

2. Experimental measurements and design studies 
of the continuous analog control MMSB system. 

3. Sampled-data analysis for digital control of 
lift and pitch axes. 

4. Hybrid simulation by using PDF 11/10 digital 
computer for control and an analog computer to simulate the 
magnetic coils and model. 

5. Description of the microcomputer. 

6. Real time digital control of the one-dimensional 
magnetic balance. 

1 . Estimation of the attainable limits of accuracy 
of model location and attitude. 

8. Conclusions and recommendations. 
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Chapter II 


DESCRIPTION OF MIT-NASA MAGNETIC MODEL 
SUSPENSION AND BALANCE SYSTEM 


Background 


MIT-NASA magnetic suspension and balance system at the 
Aerophysics Laboratory has five controllable degrees of 
freedom for ferromagnetic wind tunnel models: lift, drag, 

slip, pitch and yaw. The model can also be rotated at 
variable speed. 

Development of magnetic suspension and balance systems 
for wind tunnels was begun at M.I.T. in 1957, shortly after 
publication of the work of Tournier and Laurenceau (3) . 
Chrisinger (1) , Tilton and Schwartz (4) and Tilton and 
Baron (5) were responsible for the initial work on the 
problem done at M.I.T. 

Work on the MIT-NASA magnetic suspension and balance 
system began in 1963 under NASA-Langley Research Center 
Contract NASA-4421. It was in late 1969 that the initial 
checkout was complete (6) . 

The present MMSB system consists of the following sub- 
systems: a magnet coil assembly, a model position sensor 

system, a compensation system and a power amplifier system. 
The magnetic coil assembly generates the necessary magnetic 
field. Model position sensors sense the position of the 
model. Compensators provide the automatic control of the 
system. Power amplifiers amplify the compensated error 
signals to drive the magnetic coil assembly. 

The position sensing device (6), Electromagnetic 
Position Sensor (EPS) , is basically a multi-component 
differential transformer with the model forming the core. 

A simplified schematic is shown in Figure 2 and the proto- 
type EPS system transducer coil assembly is shown in 
Figure 3. The excitation windings are powered by a constant 
voltage amplifier at 20 KHz. Movement of the model varies 
the coupling’ between the excitation windings and the pickup 
windings. The pickup coil voltages are then processed, - 
including filtering, amplifying, summing and demodulating, 
to give high level outputs which are approximately propor- 
tional to the displacements of the model relative to the 
axes of the transducer assembly. 
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A functional block diagram of the MIT-NASA MMSB system 
is presented in Figure 1. Several assumptions were used in 
the design of the present continuous compensation system (6) . 

1. Each degree of freedom is uncoupled from the others. 
That is to say, we can treat them as uncoupled channels and 
design and analyze them independently, 

2. The transfer function from power amplifier input 
signal to the magnetic coil current is characterized by a 
first order lag. 

3. Magnetic force is assxmied to be proportional to and 
in phase with the magnet current. 

4. The measured model position is assumed to be pro- 
portional to and in phase with the actual model position. 


The first assumption is approximately valid for the MIT-NASA 
MMSB for small angles of attack, and is assumed again in the 
preliminary design of the sampled-data control system. All 
other assumptions are made again for the purposes of this 
chapter and will be discussed at its end. 

A good compensator should accomplish three things; 

1. Stabilization of the loop. 

2. Minimization of the integral static position error 
for given force disturbance inputs. 

3. Minimization of the integral static position error 
for given position command inputs. 

Consider the translational components (lift, drag and slip) 
first. They are essentially similar, characterized by the 
same form of system equations, with different coefficients. 
One channel is shown in simplified form in Figure 4. In 
order for the compensator to be quite simple and effective 
for a wide range of models and angles of attack, the form 
in Figure 5 is adopted. It consists of three paths: propor- 

tional, integral and third-order lead-lag network. The 
integrator is included to maintain zero position error at 
steady force input. 

The servo-loop diagram for this configuration is shown 
in Figure 6. In the plot there are five circled nimibers. 

They indicate the points at which data were taken. (This 
procedure will be explained in the next chapter.) With 
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appropriate combinations of gains and settings of poles and 
zeros of the compensator, satisfactory operation is achieved 
for a wide range of conditions. Root loci of lift and drag 
axes for a copper plated ellipsoidal core under certain com- 
binations are given in the next chapter. 

The pitch and yaw axes are slightly- different from the 
translational components, the main difference being the angular 
response to moment inputs. The model thus appears held by a 
torsional spring having natural frequencies determined by the 
magnetic and inertial properties of the model (2) . This has 
the effect of separating the model dynamic poles symmetrically 
from the origin, as will be shown in the root locus diagram 
of the next chapter. 

Another point is that the pitch and yaw position signals 
are supplemented by two laser angle-sensing elements. These 
elements supply the integral loops of the control systems 
holding the model at a set angular position. The laser in- 
tegrators were not used during the tests described below. 

Measured Performance of 
the Original Analogue System 

In order to have confidence in the analytical design of 
the digital system, the static and dynamic performance 
characteristics of the continuous analogue components were 
measured experimentally. The experiment included: 

1. Static calibration of the Electromagnetic Position 
Sensor system for lift, drag, pitch and yaw degrees of free- 
dom. 


2, Dynamic performance of the lift, drag and pitch 
degrees of freedom. 


The model used was the copper plated ellipsoidal core, as 
shown in Figure 7. Static position outputs of the EPS were 
calibrated with the aid of a set of transits. All static 
data were fed into a Hewlett Packard 240 1C integrating digital 
voltmeter and recorded on a HP J66562A digital recorder. 
Dynamic data were recorded on a Sanborn 4500 series recorder 
and were monitored on oscilloscopes. 

In this chapter these experiments and the data reduction 
procedure will be described, - Analytic models of the MMS^ 
system will then be presented and compared to the experimental 
data. 
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A, Static Calibration of the 
Electroftiaqnetic Position Sensor System 

Calibrations of the EPS output versus model position 
were made for the lift, drag, pitch and yaw degrees of free- 
dom. -Relations between output and position are linear, and 
there is some coupling among axes, as will be described below. 

The model was suspended at the center of the tunnel and 
then repositioned in one direction only (i.e., pure lift, 
pure drag and so forth) at several different distances. The 
relative displacement to the center of the tunnel was measured 
with the aid of a set of transits. Output of this axis as 
well as all other axes were recorded. Using the mathematical 
linear regression method, we have equations in the form 
y = ax + b; i.e,, EPS output in millivolts (y) versus displace- 
ment (x) . The results of the calibration could be summarized 
by a single matrix equation, y = ax + b, as given in Equa- 
tion 2.1. 
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The diagonal elements of the matrix A are static sensitivities 
of the EPS caused by model movements. All off-diagonal elements 
of the matrix are caused by the coupling between axes, and are 
much smaller than the corresponding diagonal elements. All 
channels also have offset voltages. 

As a matter of test procedure,' -the diagonal values have 
been checked each time the suspension has been subsequently 
used, and the values have remained within five per cent of the 
listed values during the entire testing period. However, the 
sensitivity is known to be somewhat dependent upon the relative 
humidity. 
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B. Dynamic Performance 

In order to obtain transfer functions of the continuous 
components of the system, the frequency response of the 
closed- loop continuous system was measured. We followed the 
assuH 5 >tion that each channel was essentially uncoupled, and 
could be treated independently. For the purpose of this 
design study we excluded the integrator as part of the com- 
pensator for lift and pitch axes, but included it for the 
drag axis, as will be explained at the end of this chapter. 

The copper plated ellipsoidal core was first suspended 
at the center of the tunnel under a certain combination of 
settings of gains and poles and zeros of the compensator. A 
sinusoidal input perturbation was then added as the position 
command signal, with frequency range from 0.04 Hz up to about 
20 Hz. Signals at five different points were measured, as 
shown in Figure 6, These signals were recorded on a Sanborn 
recorder . 

By measuring the voltages at these points the individual 
Bode magnitude and phase plots of each component were obtained. 
These data were then approximated by transfer functions to the 
extent possible. The resulting block diagrams are given in 
Figures 8, 12 and 17 for the lift, pitch and drag degrees of 
freedom. In these diagrams the transfer function of each com- 
ponent is normalized to its corresponding D.C. gain and the 
overall feed forward and feedback gains are put in separate 
blocks for convenience. The open loop frequency response is 
given in Figures 9, 13 and 18, For these plots the experi- 
mental data points are shown by dots, and the solid curves are 
the analytical results obtained using the assumed transfer 
function. The low frequency data permitted us to estimate the 
open-loop static gain. 

In addition to the frequency response correlation for each 
component the transfer function data was used to predict the 
closed- loop characteristics. The root locus diagrams were ob- 
tained, as shown in Figures 10, 14, 15 and 19. For the 
estimated open-loop static gains the closed- loop Bode plots 
were obtained, as shown in Figures 11, 16 and 20 for the lift, 
pitch and drag degrees of freedom, respectively. 

The correlation between theory and experiment is adequate 
for system design. Several comments are in order: 

1. As mentioned at the beginning of the Dynamic Performance 
Test, these experiments were done for a certain combination of 
settings of gains and poles and zeros of the compensator for' the 
copper plated ellipsoidal core. This means only that this com- 
pensation was workable; i.e., the model could be suspended with 
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reasonable characteristics. ‘ It does not mean the combination 
was the best. If we had changed the combination or changed the 
model, the data would have been changed. However, the closed- 
loop tests verified that the analytical modeling of the analog 
components was sufficiently adequate, and that was the point 
of this design study. 

2. The integrator is included in the compensator in order 
to maintain zero position error at steady disturbance input, 
and it is in parallel ‘■with two other paths. If the integrator 
gain Ki is comparable with 1, the effect of the integrator is 
only present at very low frequency. Because of this the inte- 
grators were excluded from the compensators for the lift and 
pitch degrees of freedom in this study to simplify data analysis 
for the other parts of the system. As we shall see in the next 
chapter, there will be no difficulty in adding an integrator to 
the sampled-data control system. 

3. These results are valid only for small displacements. 

The reason is that the suspension characteristics vary with 
model position, resulting in a non-linear system. This in- 
dicates that many things can be done in the future in the non- 
linear control area to increase the range of operation. 

4. As mentioned above several assimptions made in the 
original design have been modified by the above experiment. 
Specifically we obtained: 

a. An explicit representation of the power amplifiers 
which are characterized by second-order lags. 

b. An explicit representation of the transfer func- 
tions from magnetic field currents to the EPS outputs. For the 
lift and drag degree of freedom these are approximated by second- 
order integrators. For the pitch axis, there are a pair of 
complex poles with one real zero. These complex poles are the 
result of the spring effect mentioned above. The reason that 
they are not on the imaginary axis is believed due to eddy current 
losses in the model (7) . 
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Chapter liEI 


DESIGN OF SAMPLED DATA CONTROL 
FOR LIFT-AND-PITCH AXES 


The first requirement of a digital control system for 
the MMSB is to provide for automatic control. Since models 
are required to perform coning motion up to 15 Hz, the’ re- 
sultant system should have a good frequency response in pitch 
and yaw axes up to 15 Hz. 

The block diagram for the lift degree of freedom is 
shown in Figure 21, In the plot the analog system includes 
the power amplifier, magnetic coil, model and Electromagnetic 
Position Sensor. The digital system includes a digital filter 
in the feedback path and an integrator in parallel with the 
proportional signal in the feed forward path. The position 
signal from the EPS is sampled at the sampling frequency and 
fed to the digital filter which is put in the feedback path 
to reduce overshoot. The position command signal is sampled 
at the same rate, and added to the output of the digital 
filter. The sum is then passed through the parallel path, 
including proportional and integral components, to maintain 
zero position error at steady disturbance input. The resultant 
signal is fed to a digital-to-analog convertor and sample-and- 
hold circuit, and then fed to the analog system. 

From Figure 11 the closed- loop cutoff frequency of the 
lift axis for the present continuous control system is about 
40 radians per second (6,4 Hz). The input frequency will not 
exceed 15 Hz. The sampling frequency selected was 33.3 Hz 
with a sampling period of 0.03 second. This was chosen such 
that there was no aliasing, but was otherwise arbitrary. 

The second step was to get a discrete representation of 
the analog system, including the zero order hold circuit. 

This was done by taking the Z transform of the analog part, 

TS 

where Z = e , T = sampling period. We then specified the 
desired performance capability in terms of the natural fre- 
quency, damping ratio and zero locations of a continuous 
second order system, which was called the MODEL. The idea 
of constructing a MODEL is that one could optimize the closed 
loop response of the lift degree of freedom so as -to approxi- 
mate that of the MODEL, The MODEL used had a pair of complex 
poles at 0.7 ± 0.3j and a single zero at the origin '"in the Z 
plane. 

The integrator gain was then* selected based on the 

desired low frequency response when the model was subject to 
disturbance input. . = 0.5 was chosen with response time 
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about five seconds with respect to disturbance input. We 
then compensated the system using a digital filter to get 
the desired performance. A second order digital filter was 
selected for the lift axis. Using root locus technique, 
one then had a first estimate of the possible coefficients 
of the filter and the corresponding static open loop sensi- 
tivity, The selected coefficients were only a first guess. 

To refine the compensation characteristics , parameter 
optimization techniques were used. The parameter optimiza- 
tion technique used here was developed for the design of 
linear automatic control systems applicable to both con- 
tinuous and digital systems (8) . In these techniques the 
Model Performance Index is used as the optimization criterion 
because of the physical insight that can be attached to it. 

The design emphasis is to start with the simplest system 
configuration that experience indicates would be practical. 
Design parameters are specified, and a digital computer pro- 
gram is used to select that set of parameter values which 
minimizes the performance index. The resulting design is 
examined and complexity, through the use of more complex 
information processing or more feedback paths, is added only 
if the performance fails to meet operational specifications. 
System performance specifications are assumed to be such 
that the desired step function time response of the system 
can be inferred. 

Using the coefficients obtained from the root locus 
method as a first guess, and running the parameter optimiza- 
tion program to select the parameters of the second order 
digital filter, the system of Figure 22 was obtained. The 
Z-plane root locus for the system is shown together with the 
step response in Figure 23. The stability boundary in the 
Z-plane is the unit circle. The system is stable up to an 

-2 

open- loop gain of 280 sec , (The gain value is the open- 
loop gain with the integral compensator path open.) The 
integrator gain ratio was 0.5. 

The design procedure for the pitch axis is similar 
except that 

1. The compensated system should have a good closed- 
loop frequency response up to 15 Hz. As the result, 33.3 Hz 
sampling frequency was not fast enough. A new sampling rate 
of 100 Hz with sampling period 0,01 second was then chosen, 

, 2. The transfer function of the model has one pair of 

complex poles and one real zero, which result from the spring 
effect as explained earlier. The closed-loop uncompensated 
system has very low damping ratio. 

This system is shown in Figure 24. The Z-plane root locus 
for the system together with the step response is shown in 
Figure 25. 
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Chapter IV 


HYBRID SIMULATION 


The preliminary designs of the lift and pitch degrees 
of freedom were simulated using a hybrid digital/analog 
computer and step responses of the hybrid systems v;ere 
measured (10) . The digital computer used was the Aeronautics 
and Astronautics Department's PDF 11/10. The analog computer 
was TR-48 of Electronic Associates, Inc. The hybrid simula- 
tion program was written in BASIC and is listed in Appendix D. 

Upon entering the program the difference equation of the 
digital filter was first calculated from the assigned complex 
poles and zeros. With a given sampling period and an inte- 
grator gain the digital computer could sample the position 
and position command signals at the sampling rate and 
calculate a compensated position error signal. The resultant 
was then fed through an analog-to-digital convertor and 
zero-order hold circuit to the analog computer. 

The program was run on the PDP 11/10 for one degree of 
freedom at a time. Since the time required for one sample 
calculation was longer than the sampling period of the de- 
signed compensator, the analog system was slowed down to 
match the PDP 11/10. A time factor of 5 was chosen for the 
lift degree of freedom, and 10 for pitch axis. (Designed 
sampling period was 0.03 second and 0.01 second, respectively.) 

For the lift degree of freedom the analog computer 
system is shown in Figure 26. With the designed digital 
filter, step responses for several combinations of open loop 
sensitivity K and integrator gain were measured, as listed 

xn Table 1. Percentage overshoot and damped frequency were 
also calculated. Estimating the dominant mode characteristics 
gave the results superimposed on the root locus plot of 
Figure 23, as shown by the heavy dots. The step response for 
K = 150, and = 0.5 is given in Figure 27. The time re- 
sponse to a step disturbance force applied to the model is 
shown in Figure 28. In these plots the time scale has been 
corrected to read actual records. The pitch axis is similar 
to the lift degree of freedom. The analog computer diagram 
is given in Figure 29. Time responses for K = 16 and 32 ar,e 
given in Figure 30. 

Two comments are pertinent about the hybrid simulation 
procedure ; 
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Table 1 


No. 

S, 

X 

RESULT OF 

LIFT AXIS HYBRID 

SIMULATION 

®01 

Overshoot 

% 

Damped 

Frequency (rad) 

1 

0.5 

125 

21.0 

13.1 

2 

0.5 

150 

15.0 

15.7 

3 

0.5 

200 

9.7 

33.1 

4 

0.5 

280 


36.1 

5 

1.0 

150 

18.7 

17.5 

6 

3.0 

150 

26.7 

18.5 

7 

0 

150 

13.3 

15.7 
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1. Speed 


Some effort was made to write the program in Assembly 
language, which was believed to be much faster than the pro- 
gram written in BASIC. It was not completed before the micro- 
computer became available. The Intel, 8080 microprocessor 
with a machine cycle speed of 500 nanoseconds could accommodate 
the program in BASIC (9) . Peripherals of the microcomputer 
system included analog-to-digital convertor, digital- to-analog 
convertor and a teletype. Software facilities included a 
loader, dumper and other basic programs which were stored in 
programmable read-only memory. This microcomputer was 
assembled from stock parts of a term project and was returned 
to stock after completion of the experiments described in this 
report. It is described in detail in the next chapter. 

2. Noise 


Most noise was generated by the analog-to-digital con- 
vertors. This could be shown by fixing the analog input 
voltage and printing out the corresponding digital signals, 
which were noisy. 
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P®®CSDfi^G PAGE BLA^{^C NOT 
Chapter V 

DESCRIPTION OF THE MICROCOMPUTER 


The prototype microcomputer was based on an Intel 8080 
microprocessor and was assembled from components on loan 
from the M.I.T, Digital Systems Laboratory.* The computer 
shown in the block diagram of Figure 31 includes three major 
parts: 

Computer body 

a. CPU Module 

1) Intel 8080 microprocessor 

2) Intel 8224 clock generator 

3) Intel 8228 system controller 

4) Auxiliary controller 

b . Memory 

1) 2K programmable read- only-memory 

2) IK random-access-memory 

c. Input/Output and console 

1) Intel 8255 programmable peripheral interface 

2) Auxiliary I/O 

3) Interrupt handling circuit 

4) Front panel 

2. Peripherals 

a. Teletype and Teletype interface 

b. Eight bit A/D convertor with four-to-one analog 

multiplexer 

c. Nine bit D/A convertor with two sample-and-hold 

circuits 

d. Clock 

e. Eight channel digital input 


Construction of this microcomputer was carried out by 
Mr. Luh as part of Course 16.60 in the Electrical 
Engineering Department under Professor Hoo-Min Toong. 
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3 . Software library 

a. Console monitor 

b. Peripheral service routines 

c. Loader and dumper 

d. Arithmetic routines 

e . Control programs 


Throughout this section three nimber systems will be used; 
hexdecimal, binary and decimal. So as not to clutter all 
numbers with subscripted bases, the following general con- 
vention will be used: 

Hexdecimal: for address locations, contents of addresses 

and operation codes for instructions. 

Binary: for describing a single binary element. 

Decimal: for all normal referencing to quantities. 


A conversion table is given in Table 2, 


2’s complement subtractions are used by the 8080 machine. 


Table 2 


EQUIVALENT NUMBERS 


Decimal 


Binary Hexdecimal 


0 

0000 

0 

1 

0001 

1 

2 

0010 

2 

3 

0011 

3 

4 

0100 

4 

5 

0101 

5 

6 

0110 

6 

7 

0111 

7 

8 

1000 

8 

9 

1001 

9 

10 

1010 

A 

11 

1011 

B 

12 

1100 

C 

13 

1101 

D 

14 

1110 

E 

15 

1111 

F 
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The basic time cycle of the microcomputer (period of each 
state) is 500 nanoseconds. 

The components shown in Figure 31 are movinted on 
22 cards which are inserted into the slots of a portable 
box. Connections among the cards are made by wires at the 
back of the slots or through cables which connect different 
cards. The structure of the computer was limited by the 
components that were available. In order to avoid lengthy 
explanation about the theory and operation of the micro- 
processor itself. Reference 9 is recommended for those who 
are interested in the detailed circuits. (See Appendix A 
for Instruction List) 

Computer Body 

The structure of the computer body follows mostly the 
design done by Jon Valvano and Randy Haagens of the Digital 
Systems Laboratory, M.I.T. The schematic is given in 
Figure 32. 

. 1. CPU Module 

CPU module consists of an Intel 8080 microprocessor, 
an 8224 clock generator, an 8228 system controller, and an 
auxiliary controller. The auxiliary controller consists 
essentially of 74174 flip-flops, and is used to latch the 
status information and generate RDYIN signal for the 8224. 

2 . Memory 

The prototype microcomputer has 2K programmable-read- 
only memory (1702) with locations from 0000 to 07FF; IK 
random-access-memory (2102) with locations from 0800 to OBFF. 
Current memory location is decoded from the 16 bit address 
bus . 

Input and output of the memory are connected directly 
to the 8 bit, tristate data bus. Memory size could be ex- 
panded to 64 k with some modification of the decoding circuit. 

3. Input/Output and Console 

This essentially consists of two digital I/O devices, 
an interrupt handling circuit, and a front panel. The com- 
puter body has two digital I/O devices; 

a. Programmable peripheral interface. This is the 

Intel 8255, It has twenty- four pins which are 
individually programmed in- three groups of eight 
bit each for input or output. The three groups 
are , named port A,B,C respectively. 

b. Auxiliary I/O. This is an eight bit bistable latch 

(two of 7475) , and is used for output only. 
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Both devices connect directly to the tristate data bus. 

The interrupt handling circuit could handle two interrupts 
which are caused by two push-buttons on the front panel to 
perform several different front panel operations. 

The front panel includes twenty LEDS, four 7-segment 
displays, twelve toggle switches (higher eight are DATA 
SWITCHES, and lower four are FUNCTION SWITCHES) , and four 
push-buttons, as shown in Figure 33. Front panel operations 
are -very much like those of a traditional PDF 11 computer. 
They include : 

Reset 

Load address 


Deposit 


Reset the whole system 

Load high byte of the address 
Load low byte of the address 

Deposit a constant to the present address 
Decrease the address by one and deposit 
a constant 

Increase the address by one and deposit a 
constant 


Examine Examine the content of the present address 

Decrease the address by one and examine 
the content 

Increase the address by one and examine 
the content 


Start Run the program 

Single step Run the program with one instruction per 

command ' 


Each of these operations is caused by setting a certain com- 
bination of the lower four function switches and then pushing 
the corresponding push-button. "Load the address", "deposit 
a constant", and "examine the content" will cause interruption 
of the computer and these operations are done by corresponding 
interrupt service routines. 

. * 

Peripherals 

Peripherals included Teletype, eight bit A/D convertor, 
nine bit D/A convertor, clock and eight channel digital in- 
put. All these devices communicate with the computer body 
through the tristate Peripheral Data. Bus, which is terminated" 
at port B of 8255. A schematic is given in Figure 34. 


Equipment also borrowed 


from Digital Systems Laboratory. 
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Each device is operated by the computer body through 
a sequence of control signals. These signals are generated 
by dedicated peripheral service routines when they are 
called and decoded from port C of 8255 and the Auxiliary I/O. 

1. Teletype and Teletype Interface- 

Teletype was used at the system developing and program 
debugging stages. After those stages programs are stored in 
the PROM and variable coefficients are stored in the RAM, and 
the whole system is operable without the teletype. Teletype 
could be connected or disconnected to the system very easily. 

2 . Eight Bit A/D Convertor with 

Four~to~one Analog Multiplexer 

It includes an eight bit D/A convertor and a successive 
approximate register (Am 2503) . An analog four-to-one multi- 
plexer is used at the input end and controlled by the service 
routine. Four different analog channels could be sampled- 
A buffer is used at the output end to make the output tristate. 
The input range is from -2.5V to +2.5V, with corresponding 
output range from 00 to PF, 

, ^ 

3. Nine Bit D/A Convertor with 

Two Sample- and- ho Id Circuits 

A buffer is used at the input end to reduce the load of 
the D/A convertor to the Peripheral Data Bus. Nine bit 
D-type flip-flops (five of 7474) are used as latches follow- 
ing the buffer. Two sample-and-hold circuits are implemented 
at the output end of the D/A convertor and are separately 
controlled by the service routine. The output range is from 
OV to 5V, with corresponding input values 0000 to OlFF. 

4. Clock and Digital Input 

The clock is a synchronous eight-bit counter (two of 74161) . 
The input pulses are generated by the 8224 clock generator, 
with frequency reduced by an adjustable factor to fit different 
applications. The output of the counter goes into one set of 
a data selector (set A) to make the output tristate. The other 
set of input of the data selector (set B) is used as digital 
input device for the system (eight channels) , which is used to 
test different operation conditions. Peripheral service 
-routines -will be discussed in the next section. 
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The software development began as soon as the computer 
body had been built. It could be divided into several stages: 
1) Console monitor; 2) Peripheral service routines; 3) Loader 
and dumper; 4) Arithmetic routines; and 5) Control programs. 


1 . Console Monitor 


Front panel operations are perfoimied by this program. 

This is the most basic program which puts the computer body 
into operation. Basically it includes the RESET routine and 
two interrupt service routines. Interrupts are caused by 
pushing either the ADDRESS or the DATA push-buttons of the 
front panel. Flow charts are given in Figures 35, 36 and 37. 

2. Peripheral Service Routines 

These routines carry out the proper operation of all the 
peripherals, including teletype, A/D convertor, D/A convertor, 
clock and digital input. They monitor the operation of 
peripherals through the first three channels of the digital 
input and give appropriate commands through Peripheral Control 
Bus. They are strongly dependent upon the hardware structure 
of the computer. These routines are called each time some 
program wants to perform certain peripheral operations. The 
main program does not need to worry about the detail operations 
of the peripherals again. 

a. TTYIN: The computer receives a character from 

the teletype and puts it into the 
accumulator. 

b. TTYOUT: The computer sends the character in the 

accumulator to the teletype. 

c. A/D: Select an analog channel, convert the 

signal”to”digital form and store in the 
accumulator. 

d. D/A: The 9 bit data, which is stored in the 

lower part of the register pair HL, is 
sent to the D/A convertor. The selected 
sample-and-hold circuit then samples and 
holds the resultant signal. . 

e. Time: Read the clock into the accimulator . 

f. RESCLK: Reset the clock. 

g. TEST; Check the selected channel of the digital 

input device to see whether it is 1 or 0. 
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3, Loader and Duitiper 

These are design and debugging tools. The loader could 
load a program from punched paper tape into the RAM with a 
specified starting address. On the other hand the dumper 
dmps -the content of the specified memories (include PROM 
and RMI), onto the teletype and paper tape. By, using loader/ 
dtiraper and front panel operations a program could be written 
and debugged very easily. The flow chart of the loader and 
the dxmper are given in Figures 38 and 39. 

4 . Arithmetic Routines 

Since it is essential to have an accuracy more than eight 
bit in the position control application and the 8080 is an 
8 bit machine, double precision arithmetic is adopted in most 
cases. Floating point opera txons are desired, but they are 
too slow when done with software; therefore, fixed point 
arithmetic is used. Several subroutines were written in 
order to perform addition, subtraction, multiplication and 
other arithmetic operations. 

a. Addition, subtraction, shift left and shift right 
(with or without sign extension) 

These operations are expanded into double precision 
arithmetic in a pretty straightforward way. 

* 

b. Unsigned multiplication: BC DE * L 

Data stored in the register L and the register 
pair DE are treated as unsigned numbers; i.e., 
the eight bit content of L represents a number 

p 

from 0 to 2 -1 (255) , and the sixteen bit content 

of DE represents a number from 0 to 2^^-l. The 
contents of DE and L are multiplied together and 
the result is stored in the register pair BC. We 

know’ 2^^ * 2^ - 2^^, the result should be of 24 bits. 
However, BC is only 16 bits and we do not want to go 
into a triple-precision mode either. As- the result, 
it adds a limitation to this multiplication program: 
the result should be less than 16 bits; otherwise, 
overflow would occur. The flow chart is given in 
Figure 40. - 


*B,C,D,E,H,L are registers of 8080 and can be used as register 
pairs as BC, DE, HL. 
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5 , Control Programs 

The control programs are exactly the same for both lift 
and pitch degrees of freedom. Integrators are not included 
for the sake of simplicity. From Figure 21 the function of 
the computer is: 

a. Reset the clock and sample the model position 

signal (y) from the A/b convertor. 

b. Calculate the result of the second order differ- 

ence equation: 

z {n)-^b^y (n)+z (n) , where 

z* (n) (n-1) +a 2 Z (n-2) +b^y (n-1) ^-b^y (n-2) 

and z (n) has been calculated in the previous 
cycle (step 7) 

c. Sample the position command signal (y) and calcu- 

late the compensated error signal: e-^r-z (n) 

d. Sent output (e) through the D/A convertor. 

e. Update data: Present state becomes previous 

state .... and so forth; that is, 

z{n-l) z(n) 

z(n-2) ^ z(n-l) 


y{n-2) -f- y(n-l) 

f. Calculate z*{u) 

z* (n)-^a^z (n-l) +a 2 Z (n-2) +b^y (n-1) tt' 2 y (^“'2) 

g. Check time and go back to step 1. 


The flow chart of the program is given in Figure 41. 
Several things are worth mentioning: 
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1) Scaling of numbers: all inputs (y and r from 

A/D convertorl are 8 bit long and the coefficients and 

are scaled to 8 bit long too. With scaling factor equal to 
16, the binary nxraiber 00000001 is equivalent to 0.0625 of the 
real nvimber. Scaled coefficients of the compensators for 
lift and pitch degree of freedom are given in Table 3. 

2) Since the input range of the A/D convertor is 
from -2.5V to +2.5V, and output range of it is from 00 to FF, 
a DC constant (80) Hex is inherently added to all the signals 
by the convertor. Instead of taking this constant out 
immediately after conversion, we do something else. We know 
that the value we really want is: 

z (n)=a^z (n-D+a^z (n-2)+bQy (n)+bj^y (n-1) +b 2 y (n-2) 


However, 

a^ (z (n-l)+ 80 )+a 2 (z (n- 2 )+ 80 )+b^ (y (n)+ 80 )+b^ (y (n-l)+ 80 ) 
+£>2 (y (n-2)+80)- (a^+a 2 +bg+b^+b 2 ~l) *80=z (n)+80 


Instead of s-ubtracting the quantity 80 right after A/D con- 
version, we could have every variable biased with this DC 
mamber. In this fashion all inputs are forced to be positive 
numbers, which will simplify the use of unsigned multiplica- 
tions, The constant 80 will be cancelled out when z(n)+80 
is subtracted from the position command, which is r+80, to 
get the position error signal e. 

The time required for each complete calculation; i.e,, 
the time from the first clock-reset to the next clock-reset, 
is about 4.4 milliseconds. Since the sampling frequency is 
33.3 Hz for the lift degree of freedom and 100 Hz for the 
pitch degree of freedom, the two control programs could be 
run simultaneously, with the arrangement shown in Figure 42. 

Hybrid Simulation in Real Time 
Using the Microcomputer 


• After completion of the minicomputer the hybrid 
simulations which had been previously conducted at reduced 
speed were repeated in real time using the microcomputer 
instead of the PDP 11/10. In order to obtain the necessary 
speed, calculations were performed in fixed point arithmetic 
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Lift Pitch 


Table 3 


APPROXIMATION FOR THE COEFFICIEHTS OF THE DIFFERENCE EQUATION 
z (11) =a^z (n-1) +^ 22 ; (n-2) +b^y (n) +b^y (n-1) +b 2 y (n-2) 


Coefficient 

Designed 

Value 

Scaled Value 
in Decimal 

Scaled Value 
in Hex 

Decimal Number which the 
Scaled Niimber Represents 

^1 

.84019 

13 

OD 

.81250 


.09856 

2 

02 

.125 

■H ^0 

1.09073 

17 

11 

1.0625 

K 

^1 

1.05239 

17 

11 

1.0625 

b„ 

.22003 

4 

04 

.25 


^1 

.22024 

4 

04 

.25 

^2 

.15984 

3 

03 

.1875 

^0 

8.0853 

129 

81 

8.0625 

^1 

8.8138 

141 

8D 

8.8125 

^2 

2.1086 

34 

22 

2.125 



and the integrator in the foiiward path was eliminated. This 
produced essentially the same step responses in real time 
that had previously been obtained at reduced frequency. In 
these real time tests the sampling frequencies were 33.3 Hz 
in lift and 100 Hz in pitch. 

Step response of the lift degree of freedom for K=150 
is given in Figure 43. The position signal is pretty noisy, 
which might be caused by the following reasons: 

1. Quantization errors of coefficients 

2. Truncation errors caused by fixed point arithmetic 

3. Noise from A/D convertor 


* 

The programs are listed in Appendix B for both the lift and 
pitch degrees of freedom. 
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Chapter VI 


A DIGITAL CONTROL DEMONSTRATION USING 
THE ONE-DIMENSIONAL BALANCE SYSTEM 


Because the MIT-NASA prototype system is very complex 
and is also in frequent use for aerodynamic tests, it was 
decided to conduct initial tests of the digital control 
system using the one-dimensional demonstration balance. 

This magnetic balance is completely independent of the 
other magnetic balance systems and was used to demonstrate 
the ability of the microcomputer to perform actual digital 
control of a magnetic balance. The control software was 
rewritten to conform to the Aerophysics Laboratory one- 
dimensional magnetic balance system which has been used as 
a demonstration unit. 

The schematics of the coil and the system are shown in 
Figure 44. In this case the magnetic field was generated 
by a single coil, with the axis vertical. The model, which 
was a ferromagnetic ball, was hung underneath the coil, the 
lateral position of the ball being maintained on the coil 
axis by the natural stability- in this plane. The vertical 
displacement on the other hand was unstable open- loop. 
Closed- loop control was thus required. The microcomputer 
was used here as a digital compensator for the vertical 
direction. On-line implementation was made and the digital 
system performed well in all respects and indeed better 
than had the previous analog system. 

1. System Description 

The transfer function of the coil was a first order lag 
with the pole at -38.98 rad/sec. The transfer function of 
the model relating the position of the ball to the current 
in the coil was derived as follows; 


F; force, position upward 

X; position, position upward measured from center 
of coil 

Xq; nominal equilibrium position 
AX; X-Xq 

I; current through the coil 
Iq; nominal operating current 


?AGB BtAm not 
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K; a proportionality constant 
M; mass of the ball 


From the experiment we know: P = — =■ 6,1 

X 


Af = 


2KIAI 

X2 


2KI^AX 


= MAX 


6.2 


2KI' 


Ax + 


Mx: 


AX = 


2KI, 


Mx: 


AI 


6.3 


In the prototype case the ball is suspended below the coil 
with the nominal current required to overcome the gravity 
force. Iq = 90 amp, X^ = -.084m, and the transfer function 


Ax 

AI 


8.0 * 10 "^ 

(1-T^^) (1+ ® 


15.26' 


15.26 


) 


(m/amp ) 


6.4 


The sign of the "spring" coefficient is negative, leading to 
the unstable mode. 

We also notice that in the MIT— NASA magnetic suspension 
and balance system the corresponding transfer function for the 
lift degree of freedom is a double integrator, l/s^. The 
reason for the difference is as follows: in the MIT-NASA 

system separate coils are provided to generate a uniform 
model magnetizing field and a separate gradient field. The 
magnetic supporting force is proportional only to the gradient 
field; i.e., F is proportional to the current through the 
gradient coil only. In the demonstration case there was only 
one coil which generated both the magnetizing and the gradient 
fields , and fprce was proportional to the square of the current 
and inversely to the square of the distance. 

In addition to the coil, the ball and the microcomputer 
several other components were included in the system: 
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Position sensor — It consisted of a light source and a 
photocell. The operation of this position sensor is 
best explained by Figure 45. The output of the 
photocell was proportional to the illtuijinated area, 
and varied froiij zero to two hundred 'millivolts. To 
obtain a ± directional reference, ;a bias signal was 
provided . .• ' 

Power amplifier — This was a Torque System Incorporated 
Model PA-601 operational power amplifier. The transfer 
function of the power amplifier was a pure gain, which 
was set to a value of two. 

Generator — This was a General Electric 2CM73B7 generator, 
running at 3500 RPM. With field voltage as the input 
and armature voltage as the output, the transfer function 
was 


G(s} 


2.86 


1 + 


8.39 


volts per volt 


6.5 


Interface box — The input range of the eight bit A/D con- 
vertor was from -2.5 volts to +2.5 volts; the output 
range of the nine bit D/A convertor was from 0 volt to 
5 volts; and the output of the photocell was from 0 volt 
to 0.2 volts. Interface circuits were needed in order 
to match the different ranges of different components. 

A DC bias voltage was also necessary to generate the DC 
bias current to hold the model at the operating point. 

The interface will be described in more detail in 
Section 4 . It should be noticed’ that these interface - 
circuits include only simple algebraic operations upon 
signals, such as addition of a constant, subtraction of 
a constant, or multiplication by a constant. There was 
no phase lead or lag. As the result, in the design of 
digital compensators we did not need to consider these 
circuits at all. The block diagram of the system is 
shown in Figure 46. 

2 . Synthesis of Compensators 

After the transfer function of each component had been 
identified, the synthesis of a digital compensator follows 
the method described in Chapter III. The compensator is ' 
put in the feedback path.. Two' third-order compensators were 
designed as listed below: 
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a) zeros at 0,76427, 0,86728, .0,87893 

poles^ at 0,13534, 0,25, 0,,3, sasplin^ frequency 
~ 100 Hz and without integrator 

b) zeros at 0.68008, 0,68008, 0.79161 

poles at -0.13671, -0,18821, -0.25209, sampling 
frequency = 50 Hz with or without integrator 

The Z-plane root locus for the first compensator is given in 
Figure 47. 

Hybrid simulations using the microcomputer and the 
analog computer were carried out in the same manner as those 
described in Chapter IV prior to converting the microcomputer 
to the one-dimensional balance. 

3. Structure of the Program 

A third-order compensator was rewritten as three first- 
order compensators, as follows: 


Z (z) 
y (z) 


K(z-b2) (z-b2) (z-b^) 
(z-a^) (z-a^) 


Z3(z) 
y (z) 


Z^Cz) ^^(z) z^(z) 

Z2 (z) z^ (z) y (z) 


(z-a^) (z-a2) 


Ki (z-bi) 

(z-a^) 


and 

(n) (n- 1) +B^y (n) +B^y (n-1) 

Z2 (n)=A2Z2 (n-l)+B2Z^(n)+B3Z^ (n-1) 
Z3 (n)=A3Z3 (n-l)+B^Z2 (n)+B^Z2 (n-1) 


6.6 


6.7 


6.8 


6.9 


In the control program the compensator was calculated as 
three first-order compensators as above, with all the 
coefficients stored in the Random Access Memory. 


38 



I,f the inte^ra,tor was added in the feed- forward path 
to reduce the long term driftin^^, as shpwn in figure 46, it 
could also be treated as. a first order difference equation 
as follows : 


f 

e 


K.T . -1 

(1 t ^ 


1-z 


-1 


)) 


iC.T 

(1 + -^) 


K.T 



-1 

z 


6.10 


- Z^{n) = f(n) = A^Z^(n-l)+Bge(n)+B.^e{n-l) 


6.11 


where : 


A. 1, B — 1 + 


K.T K.T 

3 - R „ _J: ^ 

„ , - —X— i 


6.12 


The flow chart of the program is given in Figure 48. 

Whether the integrator was used or not was set by the 
second DATA SWITCH on the front panel. It was monitored 
by the program continuously. 

The control program is listed in Appendix C. 

i 

4. Interface Circuits' Control 

As described in Section VI-1, three interface circuits 
were used in order to match different operating ranges and 
to add the bias current to the coil. The open loop gain was 
also adjusted in one of the circuits. Schematic of the 
system including interface circuits is given in Figure 49. 

a) Interface A (IFA) : this received signal y' from the 

position sensor and gave the resultant signal y to the A/D 
convertor . 

Input range, (y'): 0 to .2 volts 

Desired output range (y) ; -2.5 volts to -4-2.5' volts 

Transfer function: y' = -25 (y - 0.1) 
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b) Interface R (jlFB.). ? this ^received the compensated 
error signal from the D/A convertor (il a,nd s.ent the re-> 
sult^nt signal to the power amplifier « It also included 

the open loop gain adjustment a,nd added a DC bias voltage to 
the signal to generate a bias current. 

Input range {£) : 0 to 5 volts 

Desired DC bias: 3.6 volts 

.Transfer function: f’ = - (K(f-2.5)+DC bias) 


c) Interface C: it is an input s\nnming stage. It 

could sim two input signals and also include an ON-OFF 
toggle switch. 

It was possible to operate the entire system success- 
fully. Suspension was stable and measurements were made 
of response to step input and the response to sinusoidal 
perturbation signals. Model position (y*) compensated 
error signal (f) as well as coil current (i) were recorded. 
Several step responses with sampling frequency equalled 
to 100 Hz are given in Figure 50. 

5* Points Suggested by the Experiment 

As a consequence of the experimental program there is 
increased confidence in the feasibility of digital controls 
with their many advantages. At the same time several things 
become apparent which are not likely to have been discovered 
with only theoretical analysis. These are summarized as 
follows : 

a) A fast processor is essential at the beginning if 
non-linear and adaptive controls are to be considered, since 
a high sampling frequency will be needed. Thus the computa- 
tion time available is short and more computation than used 
in this investigation will be needed. 

b) A floating point processor is desirable if the con- 
trol range is to be increased and if non-linear and/or 
adaptive controls are needed. This is necessary to avoid 
scaling limitations. 

c) In the design of digital compensators phase lead 
should be minimized to reduce the amplification of noise 
introduced in the measurement and sampling operations. 

d) Means of recording, controlling and documenting 
software are essential. 
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Chapter VII 


ATTAINABLE LIMITS OF ACCURACY 
OF MODEL LOCATION AND ATTITUDE 


The attainable limits of accuracy of -model location 
and attitude depend upon the limits of analog-to-digital 
convertors, digital~to-analog convertors, the way of 
arithmetic calculations and the word length of the 
processor used. In this chapter a rough estimation is 
made by assuming a typical minicomputer with 16 bit word 
length is used. 

From the calibration data of Equation 2.1, 
sensitivities of the EPS are: 

lift 337.2 mv/inch 

drag 540.0 mv/inch* 

pitch 16.4 mv/degree 

yaw 27.8 mv/degree 

2000 mv of EPS output, which allows movement of 5,93 
inches in lift axis, 3.70 inches in drag axis, 122.0 de- 
grees in pitch and 71.9 degrees in yaw is temporarily 
assumed to be good enough for the present analysis. 

If an eight bit analog- to-digital convertor were 
used, with the full range from zero to 2000 mv, one bit 
would represent '7.8125 mv, or 2.3169 x 10“2 inches for 
lift degree of freedom. Table 4 is' thus obtained by 
similar reasoning. With 12 bit analog— to-digital con- 
vertors, which are commercially available at reasonable 
prices, one could have fairly good resolution. 

Either ' fixed point arithmetic or floating, point 
operation could be used. With a steady state sensitivity 
of the digital filter equal to one, the digital processor 
xs assiimed to be able to calculate the compensated error 
signal without overflow. We also assume the truncation 
error caused by the finite word length is not serious. 

Assuming a 12 bit analog-to-digital convertor is 
used, with 0,488 mv per bit, the maximum quantization 
-error is 0.488 mv. To get a brief idea about the propaga- 
tion of quantization error, please refer to Figure 51. 
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‘Table 4 


RESOLUTION FOR DIFFERENT A/D CONVERTORS 


to 


Number of 


: of A/D 

mv per 

Bi-t 

Lif± 

(±n) 

8 

7.. 81 


2.32 

o 

1 

to 

10 

X.95 


5.79 

X 10“^ 

12 

4.88 X 

10'“^ 

L.45 

1 

o 

16 

XO'5 X 

X0~^ 

9.. 05 

X 10"''^ 


Corresponding Movement 


Drag 

(in) 

Pitch 

(deg) 

Yaw 

(deg) 

1.4'5 

X 10"^ 

4.67 

X 

lO"^ 

2.81 

X lO”^ 

3.61 

X XO**^ 

1.19 

X 

10“^ 

7.03 

X 10“^ 

9,. 03 

X 10”'^ 

2,. 78 

X 

10“^ 

1.76 

(-> 

o 

to 

5.6‘5 

X lO"^ 

X.B6 

X 

10“^ 

i.ro 

:x XO"^^ 



In Part a the quantization error is treated as an independent 
source, added to the system at the analog- to-digital convertor. 
The system could also be plotted, as shown in Part b. This 
system has identical eigen-values to the one ,in Part a and 
the only difference is the position of closed loop zeros. 

,As a result the responses of the two systems, are approximately 
the same. The maximum, quantization is 0.488 mv; the position 
error caused by this error is of tlie same level and is a very - 
small quantity when translated to inches or- degrees. 

It thus appears that the positional accuracy of a 
digitally-controlled MI^SB will be limited by the analogue 
position signals generated by the EPS system if the computer 
is correctly designed. 
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Chapter PAGS BLAKg 


CONCLUSIONS AND SgCOMMENDATIONS 


The resiil'ts of this feasibility .study indicate that 
present day digital computers can ' successfully perform the 
, magnetic balance control functions-. 

Actual operation of a digitally controlled one degree 
of freedom magnetic suspension system was demonstrated using 
a microcomputer constructed around the Intel 8080 processor, 
^ile sufficient for this feasibility study this processor 
is much too small for control of the full NASA prototype 
system. Even more capacity will also be required in order 
to perform data reduction and adaptive or non-linear control 
functions. 

The initial experiments described have clearly 
demonstrated the utility of a digital microprocessor for 
performing the^ magnetic balance compensation functions in 
a linear control system. They have also indicated that use 
of a digital control system should lend itself to; 

Automatic data reduction. 

Accurately preprogrammed model motion. 

Non-linear and self-adaptive control processes 
which might greatly improve performance 

The next logical step is to incorporate a dedicated digital 
processor controller into the NASA prototype magnetic balance 
system and use it to develop the above areas. After software 
has been developed to handle the desired control loop and 
data reduction functions, a system could be specified that 
would be suitable for the 2.5m NTF or pilot . 3m facility. 

One possible configuration is shown in Figure 52. The digital 
processor is divided into two subsystems, such that each of 
them could be installed using current computers at reasonable 
prices and the relationship between them is clearly defined. 

A dedicated digital processor is, used to do all the 
automatic control functions, including control of six channels 
magnetic model suspension and balance system and the manage- 
ment of the wind tunnel. The other computer is used to 
collect and reduce experimental data and carry out research 
calculations involving aerodynamic experiments and the MMSB 
itself. The second computer should be operated at high level 
language for ease of communication with tunnel operators. It 
should also have adequate memory for. future flexibility. 
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Another suggested structure is shown in Figure 53, 

One processor is used to do all the model control functions j 
including non-linear and adaptive control. The other pro- 
cessor Is used to collect and reduce experimental data and ■ 
manage the wind tunnel housekeeping. 

Another option is that the data could be stored on 
disks and analyzed later, and only one processor is needed. 
A checklist for the selection of computers is given in 
Table 5. 
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Table 5 


CHECKLIST FOR 


CPU 


Peripherals 


Software 


Future Expansion 


Service 


THE SELECTION OF COMPUTERS 


1. Architecture 

2 . Speed 

3 . Word Length 

4 . Instruction Set 

5 . Multi-Level Interrupt 

6. Direct Memory Access 

7. Floating Point Capability 

8. Reliability and Maintainability 

9 . Memory Management 
10. (Cache Memory) 

1. Terminal a, Teletype/DEC Writer 

b. CRT 

2. Disk or Fixed Head Disk 

3. Real Time Clock 

4. A/D, D/A 

5. (Line Printer) 

6 . (Floppy Disk) 

1. Real Time Operating System 

2. Language Processors 

3. Diagnostic Software 

4. (Multi-task, Multi-programming 

Operating System) 

1. Data Reduction 

2 . Tunnel Management 

3. Multi-processor System 


Price 
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Functional block diagram of suspension system 
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Figure 

Servo-loop diagram of suspension system 
translational axis (lift, drag, slip) 



Figure 

Schematic of compensator 
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Force Disturbance Input 



Figure 6 

Servo loop diagram of suspension system, translational 
degree of freedom 
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Lift degree of freedom closed-loop Bode plot 
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Figure 20 

Drag degree of freedom closed- loop Bode lot, 2/1 
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Figure 22 Lift degree of freedom sampled data analysis 
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Figure 23 
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Figure 24 Pitch axis sampled data analysis 








Normalized response 



Figure 25 Pitch axis z-plane root locus and step response 
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Figure 26 

Lift degree of freedom hybrid simulation, 
analog computer system, time factor = 5 









Figure 29 

Pitch axis hybrid simulation, analog system, 
time factor = 10 
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Figure 30b Hybrid simulation of the pitch degree of freedom. 
Step response, K=32 
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Figure 31. Block diagram of the microcomputer 
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Figure 33 Schematic of Front Panel 




Figure 34 Schematic of peripherals 
















Figure 35 

CONSOLE MONITOR (I) RESET 



Figure 36 

CONSOLE MONITOR (II) INTERRUPT (ADDRESS) 
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Figure 37 

CONSOLE MONITOR (III) INTERRUPT (DATA) 


84 









Figure 38 
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Figure 39 
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Figure 41 Plow chart of the control program 
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Figure 42 Scheduling of the microcomputer 
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Figure 46 Block diagram.' of the one-dimensiohal Balance system 










Figure 47 Z~plane root locus for 100 Hz sampling frequency 
digital compensator for one-dimensional balance 



Pigure 48 Flow chart of the program for the 
one-dimensional balance system 
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Figure 49 Schematic of the one-dimensional balance system, 
including interface circuits 










ORIGINAL PAGE IS 
OE POOR QUALITY 


Model 



Qaantizatdon error 


(a) 


Position cxnmand 











VO 

03 



Figure 52 

Suggested structure of 
the digital processor, I 










Figure 53 

Suggested structure of the 
digital processor, II 







Appendix A. INSTRUCTION LIST FOR INTEL 8080 MICROCOMPUTER 


»ilUMP 


CALL 


RETURN 


RESTART 


ROTATSt 


MOVE (coni) 

ACCUMULATOR* 




CONSTANT 






















AS 

XRA 

8 

DEFINITION 

C3 

JMP 


CD 

CALL" 


C9 

RET 


C7 

RST 

0 

07 

RLC 


58 

MOV 

E.B 

80 

ADD 

B ' 

A9 

XRA 

C 



C2 

ONZ 


C4 

CNZ 


CO 

RNZ 


CF 

RST 

1 

OF 

RRC 


59 

MOV 

E.C 

81 

ADD 

C 

AA 

XRA 

0 

obdhI 

1 

xCA 

JZ 

1 

CC 

CZ 


C8 

RZ 


D7 

RST 

2 

17 

HAL 


5A 

MOV 

e.D 

82 

ADD 

0 

AB 

XRA 

E 

1AHj 

* hw 

1 

02 

Jnc ’ 

1 

1 

04 

CNC 


DO 

RNC 


OF 

RST 

3 

IF 

RAH 


58 

MOV 

E.E 

83 

ADO 

E 

AC 

XRA 

H 



IDA 

JC 

Vaw, DC 

CC 

> Adr 

08 

RC 


» E7 

RST 

4 




5C 

MOV 

E,H 

84 

ADD 

H 

AD 

XRA 

L 

1050 1 

' f-k 

E2 

JPO 

r 

E4 

CPO 


EO 

RPO 


EF 

RST 

5 




5D 

MOV 

E,L 

85 

ADD 

L 

AE 

XRA 

M 

105 

> DQcirnai 

£A 

JPE 


EC 

CPE 


'E8 

RPE 


F7 

RST 

e 




5E 

MOV 

E.M 

86 

ADD 

M 

AF 

XRA 

A 

1 



JP 


F4 

CP 


FO 

RP 


FF 

RST 

7 

CONTROL 


5F 

MOV 

E.A 

67 

ADO 

A 




720 1 


■"FA 

JM , 


FC 

CM ^ 


F8 

RM 














BO 

ORA 

B 

720 

> ucidi 

fE9 

PCHL 











00 

NOP 


60 

MOV 

H.B 

68 

ADC 

B 

B1 

ORA 

C 


1 













76 

HLT 


61 

MOV 

H.C 

89 

ADC 

C 

B2 

ORA 

D 

iioiib1 














F3 

Dl 


62 

MOV 

H.D 

8A 

ADC 

D 

B3 

ORA 

E 

0011 OB 

> B/na?y 













FB 

El 


63 

MOV 

H,E 

8B 

ADC 

E 

B4 

ORA 

H 


J 

.MOVE 


•Acc 



LOAD 







• 

64 

MOV 

H,H 

80 

ADC 

H 

B5 

ORA 

L 

TEST' 


StMMEOIATE 


IMMEDIATE’ 

IMMEDIATE 


STACK OPS 




65 

MOV 

H.L 

8D 

ADC 

L 

B6 

ORA 

M 

•A’ ‘8’ 

ASCH 
















66 

MOV 

H.M 

8B 

ADC 

M 

87 

ORA 

A 



*06 

MVI 

B, 'I 

C6 

ADI " 


01 

Ui 

B. 

05 

PUSH 8 



67 

MOV 

H,A 

8F 

ADC 

A 






IDE 

MVI 

c, • 

ce 

AC! 


11 

LXI 

D. 

. D5 

PUSH D 










B8 

CMP 

B 

OPERATORS 

’-16 

MVI 

D. 

D6 

SUI 


21 

LXl 

H, 

ID16 £5 

PUSH H 

40 

MOV 

8,B 

68 

MOV 

L.B 

90 

SUB 

3 

B9 

CMP 

C 



'.t£ 

MVI 

E. 

DE 

S8i \ 


31 

Uf 

SP. 

F5 

PUSH PSW 

41 

MOV 

B.C 

, 69 

MOV 

L.C 

91 

SUB 

C 

BA 

CMP 

D 



.26 

MV! 

H, 

E6 

ANI 

► 08 







42 

MOV 

B.D 

6A 

MOV 

UD 

92 

SUB 

D 

BB 

CMP 

E 

• / . MOD.SHL.SH(- 

"2E 

MVI 

L, 

EE 

XRI 





Cf 

POP 

B 

43 

MOV 

B,E 

6B 

MOV 

L,E 

93 

SUB 

E 

BC 

CMP 

H 




MVI 

M, 

F6 

ORI 





D1 

POP 

D 

44 

MOV 

B,H 

6C 

MOV 

L.H 

94 

SUB 

H 

BO 

CMP 

L 

NOT 

o o 

:S£ 

MVt 

A. J 

FE 

CPI ,/ 





El 

POP 

H 

45 

MOV 

B,L 

60 

MOV 

L,L 

. 95 

sue 

L 

BE 

CMP 

M 

AND 








-DOUBLE AOD+ 

FI 

POP 

PSW' 

46 

MOV 

B,M 

6E 

MOV 

L,M 

96 

SUB 

M 

8F 

CMP 

A 

OR.XOR 


:2> 






09 

DAD 

B 




47 

MOV 

B.A 

6F 

MOV 

UA 

97 

SUB 

A 












19 

DAD 

0 

E3 

XTHL 


48 

MOV 

C,B 

70 

MOV 

M,B 

98 

SBB 

B 












29 

DAD 

H 

F9 

SPHL 


49 

MOV 

C C 

71 

MOV 

M.C 

99 

SBB 

C 

PSEUDO 


STANDARD ® ^ 

iINCREMENT” 

DECREMENT”- 

39 

DAD 

SP 




4A 

MOV 

N/jW 

c.o 

72 

MOV 

M.D 

9A 

SBB 

0 

fNSTRUCTfON 

SETS 














46 

MOV 

C.E 

73 

MOV 

M,E 

SB 

SBB 

£ 





F S 

i04 

INR 

B 

05 

OCR 

B 

\ 



SPECIALS 


4C 

MOV 

C.H 

74 

MOV 

M.H 

90 

SBB 

H 

ORG Adr 


A SET 7a 

yoc 

INfl 

C 

OD 

OCR 

C 







4D 

MOV 

C.L 

75 

MOV 

M.L 

90 

SBB 

L 

END 


B SET 0r< ^ 

T>4 

(NO 

D 

15 

OCR 

D 

LOAO/STORE 

EB 

XCHG 


4E 

MOV 

C.M 




9E 

SBB 

M 

EQU D16 


C SET '-f 

iC 

INIR 

E 

ID 

DCR 

E 




27 

DAA* 


4F 

MOV 

C.A 

77 

MOV 

M,A 

9F 

S8B 

A 

SET 

D16 


D SET 2 

:24 

INR 

H 

25 

OCR 

H 

OA 

■LDAX 

8 

2F 

CMA 











DS 

D16 


E SET 3 

2C 

INR 

L 

'2D 

DCR 

L 

1A 

LDAX 

D 

37 

STCt 


50 

MOV 

D.a 

78 

MOV 

A.B 

AO 

ANA 

B 

DQ 

08 

11 

H SET 4 

134 

INR 

M 

35 

DCR 

M 

2A 

LHLD Adr 

3F 

CMCt 


51 

MOV 

o.c 

79 

MOV 

A.C 

A1 

ANA 

C 

OW 

'016 

( 1 

L SET 5 

,3C 

INR 

A 

3D 

DCR 

A 

3A 

LDA 

Adr 




62 

MOV 

D,D 

7A 

MOV 

A.D 

A2 

ANA 

D 




M SET 6 













53 

MOV 

D.E 

7B 

MOV 

A.E 

A3 

ANA 

E 

IF 

D16 


SP SET 6 

03 

(NX 

B 

08 

OCX 

B . 

02 

STAX 

B 

INPUT/OUTPUT 

54 

MOV 

D,H 

7C 

MOV 

A,H 

A4 

ANA 

H 

ENDIF 


PSW SET 6 

13 

INX 

0 

18 

OCX 

0 

12 

STAX 

0 


1 


55 

MOV 

0,1 

70 

MOV 

A,L 

A5 

ANA 

L 






23 

INX 

H 

2B 

DCX 

H 

22 

SHLD Adr 

D3 

OUT 1 

'D0 

56 

MOV 

D.M 

7E 

MOV 

A.M 

A6 

ANA 

M 

MACRO 

I ) 

FLAG BYTE 

33 

INX 

SP 

3B 

OCX 

SP 

32 

STA 

Adr 

DB 

IN 1 

57 

MOV 

D.A 

7F 

MOV 

A.A 

A7 

ANA 

A 

ENDM 


STACK FORMAT 

dlT 











j 














7 6 573210 

^0-; 

» oonalant. or logical/anihmetio oxprsaslon that evaluatas 

D16 

* constant, or logicaVarlthmetie exprasslon that evaluates 

Adr » 16 bit address 




s Vc 


to an 8 bit data guamity. 






to a 16 bit data quantity. 





” <« all Rags except CARRY affected; 



air Flags (C, Z, S, P, AC) atlacted 


t « only CARRY aWected 


(exceplion; INX & OCX affect no Flags) 



Appendix B. 


ORIGINAL PAGE IS 
OF POOR QUALITY 
PJIOGRAMS FOR LIFT AND PITCH 


Address 


OP Code Mnemonic Code 


04 

00’ 

3E 

MVI,A 


1 

FF 

• 


2 

32 

STA 

- 

3 


DD 


04 


OB 


5 

D3 

OUT 4 


6 

04 



7 

CD 

CALL I.C. 


08 


60 


9 


06 


A 

3E 

mi, A 


B 

08 



OC 

CD 

CALL TEST 


D 


30 


E 


03 


F 

C2 

JNZ* 


10 


OA 


1 


04 

ENTl 

2 

CD 

CALL RESC: 


3 


BO 


14 


03 

ENT2 

15 

3E 

MVI,A 

- 

6 

FD 



7 

CD 

CALL A/D 


18 


00 


9 


03 


A 

32 

STA 


B 


^7 


1C 


OB 


D 

5F 

MOV E.A. 


E 

16 

MVI,D 


P 

00 


04 

20 

2E 

MVI,L 


1 

81 



• 2 

CD 

CALL MUL 


3 


90 


24 


03 


5 

2A 

LHLD 


6 


D1 


7 


OB 


Commeni: 

Start 


Set up initial 
conditions 


Begin to operate 


If no, jump to * 


If yes, reset clock 


Sample y 


Z-^Z*+b y 
o 
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B-2 


Address 

04 28 

9 
A 
B 

2C 

D 

E 

P 

30 

1 

2 

3 

34 

5 

6 
7 

38 

9 

A 

B 

3C 

D 

E 

F 

40 

1 

2 

3 

44 

5 

6 
7 

48 

9 

A 

B 

4C 

D 

E 

F 


OP Code 

Mnemonic Code 

CD 

CALL ADDl 

OE 

58 

03 

MVI,C 

00 

06 

MVI,B 

07 

CD 

CALL SUBl 

CD 

60 

03 

CALLv4 

CD 

90 

06 

CALLt4 

22 

90 

06 

SHLD 

3E 

D1 

OB 

MVI,A 

PC 

CD 

CALL A/D 

32 

00 

03 

STA 

4F 

DA 

OB 

MOV C.A. 

06 

MVI,B 

00 

CD 

CALL SUB2 

22 

F8 

03 

SHLD 

CD 

DB 

OB 

CALL-S-4 

CD 

90 

06 

CALLv4 


Comment 


Sample r 


e-*-r-z (n) 
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B-3 


Address 


04 


3S 

OP Code 

Mnemonic Code 

50 


90 

1 ' 


06 

2 

OE 

MVI,C 

3 

FF 


54 

06 

MVI,B 

5 

00 


6 

CD 

CALL ADDl 

7 


58 

58 


03 

9 

3E 

MVI,A 

A 

DF 


B 

CD 

CALL D/A 

5C 


CO 

D 


02 

E 

3E 

MVI,A 

F 

10 


60 

CD 

CALL TEST 

1 


30 

2 


03 

3 

2E 

MVI L 

64 


D7 

5 

26 

mi H 

6 


OB 

7 

CC 

- CZ PRINT 

68 


CO 

9 


06 

A 

2E 

MVI L 

B 

D8 


6C 

26 

MVI H 

D 

OB 


E 

CD 

CALL MOVE 

F 


DO 

70 


03 

1 

CD 

CALL S a.Z. 

2 


AO ^ ^ 

3 


04 

74' 

3A 

LDA 

5 


- DD 

6 


OB 

7 

E6 

ANI 


Comment 


Output e 


Update data 


Calculate Z* 
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B-4 


Address 

OP Code 

Mnemonic Code 

Comment 

04 78 

FF ' 



9 

C8 

RZ 


A 

3E 

JSIVI A 

Hold? 

B 

08 



1C 

CD 

CALL TEST 

■ 

D 


30 


E 


03 

If hold, j\3mp ** 

P 

C2 

JNZ** 


80 


7A 


1 


04 


*** 2 

CD 

CALL TIME 

Time for next 




sampling 

3 


CO 


84 


03 


5 

DE 

SBI 


6 

EA 



7 

FA 

JM*** 

No jump *** 

88 


82 


9 


04 


A 

C3 

JMP ENTl 

Yes, jump ENTl 

B 


12 


8C 


04 



D 

E 

F 


90 

1 

2 

3 

94 

5 

6 
7 

98 

9 

A 

B 

9C 

D 

E 

P 


J 

i 


104 



B-5 


Address 


OP Code Mnelnonic Code 
Calculate Z* 


04 


AO 

3A 

LDA 

1 


D9 

2 


OB 

3 

5F 

MOV E.A. 

A4 

16 

IWI D 

5 

00 


6 

2E 

MVI I* 

7 

22 


A8 

CD 

CALL MUL 

9 


90- 

A 


03 

B 

60 

MOV H.B. 

AC 

69 

MOV L C 

D 

22 

SHLD 

E 


D1 

F 


OB 

BO 

2A 

LHLD 

1 


D3 

2 


OB 

3 

5D 

MOV E.L. 

B4 

54 

MOV D.H. 

5 

2E 

MVI.L 

6 

04 


7 

CD 

CALL MUL 

B8 


90 

9 


03 

A 

2A 

LHLD 

B 


D1 

BC 


OB 

D 

CD 

CALL SUBl 

E 


60 

F 


03 

CO 

. 22 

SHLD 

1 


Dl 

2 


OB 

3 

2A 

LHLD 

C4 


D5 

5 


OB 

6 

5D 

MOV E.L. 

7 

54 

MOV D.H. 


Comment 




1 X 
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B-6 


Address 

OP Code 

Mnemonic Code 

04 C8 

2E 

MVI.L 

9 

03 

a. 

A 

CD 

CALL MUL 

B 


90 

CC 


03 

D 

2A 

LHLD 

E 


D1 

F 


OB 

DO 

CD 

CALL SUBl 

1 


60 

2 


03 

3 

22 

SHLD 

D4 


Dl 

5 


OB 

6 

3A 

LDA 

7 


D8 

D8 


OB 

9 

5F 

MOVE. A 

A 

16 

lyiVL.D 

B 

00 


DC 

2E 

MVI.L 

D 

8D 


E 

CD 

CALL MUL 

F 


90 

EO 

03 


1 

2A 

LHLD 

2 


Dl 

3 


OB 

E4 

CD 

CALL SUBl 

5 


60 

6 


03 

7 

22 

SHLD 

E8 


Dl 

9 


OB 

A 

C9 

RET 


B 


EC 

D 

E 

F 


Comment 


Z* 
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B-7 


Address OP Code Mnemonic Code 


04 FO 
1 
2 
3 

F4 

5 

6 
7 

F8 

9 

A 

B 

FC 

D 

E 

F 


05 


ENTl 

ENT2 


Pitch 


00 

3E 

MVI.A 

1 

FF 


2 

32 

STA 

3 


DD 

04 


OB 

5 

D3 

OUT 4 

6 

04 


7 

CD 

CALL I.C. 

08 


60 

9 


06 

A 

3E 

MVI.A 

B 

08 


OC 

CD 

CALL TEST 

D 


30 

E 


03 

F 

C2 

JNZ* 

10 


OA 

1 - 


05 

2 

CD 

CALL RESCLK 

3 


BO 

14 


03 

15 

3E 

MVI.A 

6 

FF 


7 

CD 

CALL A/D 


Comment 


Start 


Set up initial 
conditions 


Begin to operate 


If no, jump * , 


Yes, reset clock 


Sample y 
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Bt-8 


Address 


OP Code Mnemonic ;Cpde 


ENT2 


05 


18 


00 

9 


03 

A 

32 

STAU. 

B 


F7 

1C 


OB 

D 

5F 

MOV E.A. 

E 

16 

MVi D 

F 

00 

- 

20 

2E 

MVI.L 

1 

11 


2 

CD 

CALL MUL 

3 


90 

24 


03' 

5 

2A 

LHLD 

6 


FI 

7 


OB 

28 

CD 

CALL ADDl 

9 


58 

A 


03 

B 

OE 

MVI.C 

2C 

00 


D 

06 

MVI.B 

E 

07 


f 

CD 

CALL SUBl 

30 


60 

1 


03 

2 

CD 

CALLv4 

3 


90 

34 


06 

5 

CD 

CALLv4 

6 


90 

7 


06 

38 

22 

SHLD 

9 


FI 

A 


OB 

B 

3E 

MVI. A 

3C 

FE 

* 

p 

CD 

CALL A/D 

E 


bo 

F 


03 


vComment 


Z'^z*+b y 


Sample r 
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Address 

B-9 

OP Code Llnetnonic Code 

Coiranent 

05 40 

32 

STA 


1 


FA 


2 


OB 


3 

4F 

MOV C.A, 

e r-Z (n) 

44 

06 

MVI.B 


5 

00 



6 

CD 

CALL SUB2 


7 


F8 


48 


03 


9 

22 

SHLD 


A 


FB 


B 


OB 


4C 

CD 

CALLt4 


D 


90 


E 


06 


F 

00 

NOP 


50 

00 

NOP 


1 

00 

NOP 


2 

OE 

MVI.C 


3 

FF 



54 

06 

MVI.B 


5 

00 



6 

CD 

CALL ADDl 


7 ' 


58 


58 


03 


9 

3E 

MVI.A 

Output e 

A 

BF 



B 

CD 

CALL D/A 


5C 


CO 


D 


02 


E 

3E 

MVI.A 


F 

10 



60 

CD 

CALL TEST 


1 


30 


2 


03 


3 

2E 

MVI.L 


64 

F7 



5 

26 

MVI.H 


6 

OB 



7 

CC 

CZ PRINT 
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B-10 


Address 

OP Code 

Mnemonic Code 

05 68 


CO 

9 


06 

A 

2E 

MVI.L 

B 

F8 


6C 

26 

MVI.H 

D 

OB 


E 

CD 

CALL MOVE 

F 


DO 

70 


03 

1 

CD 

CALL Z a . Z . 

1 1 

2 


AO 

3 


05 

74 

3A 

LDA 

5 


DD 

6 


OB 

7 

E6 

ANI 

78 

FF 


9 

C8 

RZ 

** A 

3E 

MVI.A 

B 

08 


7C 

CD 

CALL TEST 

D 


30 

E 


03 

C 

C2 

JNZ** 

80 


7A 

1 


05 

*** 2 

CD 

CALL TIME 

3 


CO 

84 


03 

5 

DE 

SBl 

6 

4E 


7, 

FA 

JM*** 

88 


82 

9 


05 

A 

C3 

‘ IMPEND 

B 


12 

8C 


05 


D 

B 


F 


Comment 


Update data 


Calculate Z* 


Hold? 


If hold, jump ** 


Time for next 
sampling 


No, jump *** 


Yes, jump 
ENTl 
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B-11 


Address OP Code Mnemonic Code 

05 90 

1 
2 
3 

94 

5 

6 
7 

98 

9 

A 

B 

9C 

D 

E 

F 


Calculate Z* 


AO 

3A 

LDA 

1 


F9 

2 


OB 

3 

5F 

MOV E.A. 

A4 

16 

MVI.D 

5 

00 


6 

2E 

MVI.L 

7 

04 


A8 

CD 

CALL MUL 

9 


90 

A 


03 

B 

60 

MOV H.B. 

AC 

69 

MOV L.C. 

D 

22 

SHLD 

E 


FI 

F 


OB 

BO -• 

2A 

LHLD 

1 


F3 

2 


OB 

3 

5D 

MOV E.L. 


Comment 
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B-12 


Address 


OP Code 


Mnemonic Code 


05 


B4 

5 

'6 


54 

2E 

OD 


MOV D.H. 
MVI.L 


7 

CD 

CALL MUL 

B8 


90 

9 


03 

A 

2A 

LHLD 

B 


PI 

BC 


OB 

D 

CD 

CALL ADDl 

E 


58 

P 


03 

CO 

22 

SHLD 

1 


PI 

2 


OB 

3 

2A 

LHLD 

C4 


F5 

5 


OB 

6 

5D 

MOV E.L. 

7 

54 

MOV D.H. 

C8 

2E 

MVI.L 

9 

02 


A 

CD 

CALL MUL 

B 


90 

CC 


03 

D 

2A 

LHLD 

E 


FI 

F 


OB 

DO 

CD 

CALL SUBl 

1 


60 

2 


03 

3 

22 

SHLD 

D4 


FI 

5 


OB 

6 

3A 

LDA 

7 

■ 

F8 

D8 


OB 

9 

5F 

MOV E.A. 

A 

16 

MVI.D 

B 

00 



Comment 




Z*^Z*-h^U^ 
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THIS PAGE LEFT BLANK INTENTIONALLY 
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B-13 


Address 

OP Code 

Mnemonic Code 

05 DC 

2E 

MVI.L 

D 

11 

^1 

E 

CP 

CALL MUL 

F 


90 

EO 


03 

1 

2A 

LHLD 

2 


PI 

3 


OB 

E4 

CD 

CALL SU31 

5 


60 

6 


03 

7 

22 

SKLD 

E8 


FI 

5 


OB 

A 

C9 

RET 


B 

EC 

D 

E 

F 

FO 

1 

2 

3 

F4 

5 

6 
7 

F8 

9 

A 

FC 

D 

B 

F 


Comment 
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Appendix C. PROGRAM FOR ONE- DIMENSIONAL BALANCE 


Address 

OP Code 

Mnemonic Code 

Comment 

07 . 00 

3E 

MVI,A 

Start 

1 

FF 



2 

D3 

OUT 4 


3 

04 



04 

CD 

CALL IC3 

Set up initial 




conditions 

5 


AO 


6 


07 


7 

3E 

MVI,A 

Begin to opera 

08 

08 



9 

CD 

CALL TEST 


A 

30 



B 

03 



OC 

C2 

JNZ* 

If no, jump * 

D 


07 


E 


07 

If yes, reset 




clock 

ENT F 

CD 

CALL RESET 


07 10 

BO 



1 

03 



2 

3E 

MVI,A 

Sample y 

3 

FD 



14 

CD 

CALL AD 


5 

00 



6 

03 



7 

32 

STA 


18 


00 


9 


09 


A 

2E 

MVI, L 

Zi^B^y-Z^* 

B 

00 



. 1C 

CD 

CALL CAL 


D 


CO 


E 


04 


F 

2E 

MVI, L 

V®2"r"2* 





20 

06 



1 

CD 

CALL CAL 


2 


CO 


3 


04 


24 

2E 

MVI, L 

"3"®4'=2“23* 

5 

OC 



6 

CD 

CALL CAL 


7 


CO 
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C-2 


Address 

OP Code 

Mnemonic Code 

07 28 


04 

9 

3E 

MVI. A 

A 

FC 


B 

CD 

CALL AD 

2C 

00 


D 

03 


E 

32 

STA 

F 


20 

30 


09 

1 

06 

MVI, B 

2 

00 


3 

4F 

MOV C.A. 

34 

2A 

LHLD 

5 


12 

6 


09 

7 

CD 

CALL SUB2 

38 

F8 


9 

03 


A 

22 

SHLD 

B 


16 

3C 


09 

D 

3E 

MVI, A 

E 

10 


F 

CD 

CALL TEST 

40 

30 


1 

03 


2 

C2 

Jtgz** 

3 


4A 

44 


07 

5 

2E 

MVI, L 

6 

16 


7 

CD 

CALL 

48 


28 

9 


08 

A 

CD 

CALL T 4 

B 

90 


4C 

06 


D 

3A 

LDA 

E 


13 

F 


08 


Comment 
Sample r 

e-^-r-Z^ (n) 


Integrator in or 
not 

If not, jump ** 

If yes, 

Z4.^B^e+Z4* 
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Address 

OP Code 

C-3 

Mnemonic Code 

Comment 

07 50 

4F 

MOV C.A 


1 

3A 

LDA 


2 


14 


3 


08 


54 

47 

MOV B.A 


5 

CD 

CALL ADDl 


6 

58 



7 

03 



58 

3E 

MVI.A 

Output f 

9 

DF 



A 

CD 

CALL DA 


B 

CO 



5C 

02 



D 

CD 

CALL RELOCATE 

Update data 

E 


80 


F 

05 



60 

2E 

MVI.L 

Calculate 2,* 

1 

03 


1 

2 

CD 

CALL SUM 


3 


20 


64 


05 


5 

2E 

MVI.L 

Calculate Z_* 

6 

09 


2 

7 

CD 

CALL SUM 


68 


20 


9 


05 


A 

2E 

MV1,L 

Calculate Z^* 

B 

OF 


J 

6C 

CD 

CALL SUM 


D 


20 


E 


05 


F 

2E 

MVI,L 

Calculate Z^* 

70 

19 



1 

CD 

CALL SUM 


2 

20 



3 

05 



*** 74 

3E 

MVI,A 

. Hold or not? 

^ 5 

08 



6 

CD 

CALL TEST 


7 

30 
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C-4 


Address 

OP Code 

Mnemonic Code 

*** 

78 

03 


JL 

9 

C2 

JNZ*** 


A 


74 


B 


07 


7C 

3A 

LDA 


D 


12 


E 


08 


P 

47 

MOV B.A. 

07+ 

80 

CD 

CALL TIME 


1 

CO 



2 

03 



3 

90 

SUB B 


84 

FA 

JM+ 


5 


80 


6 


07 


7 

C3 

JMP ENT 


88 


OF 


9 


07 


A 

B 


8C 

D 

E 

F 

07 90 

1 
2 
3 

94 

5 

6 
7 

98 

9 

A 

B 

9C 

D 

E 

P 

0 

1 

2 

3 


Comment 


If yes, jump 


Check Time 


Jump ENT 


*** 
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C-5 


Address 

OP Code 

Mnemonic Code 

Comment 


07 4 

5 

6 

7 

8 
9 
A 
B 

C 

D 

E 

F 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 
A 
B 

C 

D 

E 

F 


08 00 
1 
2 
3 

04 

5 

6 
7 

08 

9 

A 

B 


Memory layout of coefficients 





00 


00 
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C-6 


Address 

08 tic 

D 

E 

F 

10 

1 

2 

3 

14 

5 

6 
7 

18 

9 

A 

B 

1C 

D 

E 

P 

20 

1 

2 

3 

24 

5 

6 
7 

28 

9 

A 

B 

2C 

D 

E 

F 

0 

1 

2 

3 


OP Code Mnemonic Code 


CD 

C9 


A 

A3 

Sl 

^3H 

TIME 

ADICj 

ADIC 


H 


A, 


:4l 

'4H 


ADDSUB 1 
2 

3 

4 

^^IH' ®CH 

IC^^ 

^2H 


IC 

IC 

IC 

IC 


3L 

3H 

4L 

4H 


CALL 

CO 

05 

RET 


Comment 


00 


00 

41 

3F 

40 

00 

00 

00 


00 


CD 

CO 

05 

C9 
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C-7 


Address OP Code Mnemonic Code 


ommenf 


08 4 

5 

6 

7 

8 
9 
A 
B 

C 

D 

E 

F 


Memory layout of variables 


09 


00 

1 

2 

3 


U(n) 

0 

U(n-l) 

0 


04 

5 

6 
7 


llSL 

yisn 

tUn)lL 

^KnH 


08 

9 

A 

B 


Z 

Z 

z 

z 


l{n-l)L 

l(n-l)H 

2SL 

2SH 


OC 

D 

E 

F 


l2{n)L 

p{n)E 

72(n01)L 

^2{n-l)H 


10 

1 

2 

3 


hsh 
psn 
p(n)L 
3 (n)H 


14 

5 

6 
7 


-^3(ri-i)L 

f3(n-l)H 

(n)H 


18 

9 

A 

B 


e 


e 

Z 

Z 


(n-l)L 

(n-l)H 

4SL 

4SH 
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C-8 


Address OP Code Mnemonic Code 


09 1C 
D 
E 
F 

20 

1 

2 

3 

24 

5 

6 
7 

28 

9 

A 

B 

2C 

D 

E 

F 

30 

1 

2 

3 

34 

5 

6 
7 

38 

9 

A 

B 

3C 

D 

E 

F 


^4(n)L 

Uin)E 

^4(n-l)L 

^4(n-l)H 

XI 


Comment i 
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Appendix D. BASIC PROGRAM ROR HYBRID SIMULATION 


HSPJ J5-IiFC~/7 BABIfVCAFS U0J“01 

iO RCM IIYBRJD SI HUI.A f I (IN PRO(3PAM L 
20 RPM JANUARY PlMjV/7 
30 PRINT "ENTFR NO. or 7PK0B“ 

40 INPUT N 
50 PRINT 
60 FDR I I I D N 

70 PRIN1 "LIST cXROy ONF nT A TrMF“ 

SO INPUT A<'4>^'"B( h)' ^ ’ 

90 NEXT T . V, f ^ 

LOO PRTNl' ‘ 
liO LET M-N 

T 20 GOBUB ;;:oo 

!30 FOR 1 = 0 TO N 
140 LET D( I) CCD 
345 NEXT I 
3 50 k:.: kl 

3 55 PRINT "ENILR N0<- OF >'01 IS" 

L60 INPUT N 
365 PRINT 
1 70 F OR J ■•= i 10 n 

180 PRINT "I rST PULEv ()NF Al A liME" 

185 INPUT A< L) .H( L) 

J90 NEXT I 
395 PR [NT 
200 GOBUB L20<i 
205 k-kl'^kj 

210 PRJNT '‘INTTR ORIGINAL OAMi 'L F pri\'jn\i., (il ( N; 1, tifjUN l-AEU'R" 

215 INPUT lyO 

220 PRINT 

225 T1~=TYS;T:60 

?35 FOR !“J TO h 

240 U( J / -0 

245 1(1 '-0 

250 NEXT I 

255 E>.3''-0 

260 f ' i ' -0 

265 TOR l~3 ID N 

2^0 A ( I ^ = C ( ; 3 

2 SO NFXr J 

290 PQR L-0 !0 iN 

BOO b(N“N-^L ‘ U( ! '5}!|\ 

BIO NEXT L 

B20 print "NF.FTi DllFFRENGL UTUAT LUN OR NOT'T" 

B25 INPUT J . - 

"B30 i-R[NT • ' ' 

535 II J=0 GO ro 440 

550 PRINT "TUI- U I ■ FI RENCI- I-KTIIA I I UN if. " 

BdO PRINl "ZCN)~ 

37 0 F OR I - I TO M 

380 PRINT A ( f ) « " 7 ( N - " 0 I ? " ) { " y 

390 NEXT I 

410 FOR I=H--N TO M~ I 


ORIGMAi; PAGfE It 
^OOB QlTALITig 
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420 PRINT B(I) S y [y " ) r y 

425 NEXT I 

430 PKINf B(H)y"U<N~"yMy“)" 

435 PR TNT 

440 PR[Nr "ENTER LNTEGRArOR BATN'’ 

445 LNPUr R 

450 PRINT 

455 C0-=K)KT/2+J 

460 CL=-K>!<r/2“ L 

465 AOUTd y ^65) 

470 Z(0;- 0 
475 F(0'»=0 

480 PRINT "READY 5 0 OPti^ATE NOW" 
485 AINP(3y|-l) 

490 ir il "2 GO lU 500 
495 00 TO 4H5 
500 riMR 
505 AINPC 1 vX I ) 

510 AINP(2yX2) 

515 A0Ur<.2yJ) 

520 U(0^-=X2 

525 2 < 0 ) - Z ( 0 ) HU H~N ) >KIK li-N 5 
530 L(0) -XL"-2(0) 

540 r(0)-^r(oj+co^E<o> 

550 X--F<<n-,67 
560 AOUT ( LyX) 

570 A0Ln(2y0) 

580 rOR !- ! iO h 
585 Zdl- L f i )-Z(h - ! 

590 IKh i H I / 

595 NEXT 1 
600 I 

605 1- < J > -r(0) 

610 2(0) --0 

615 I DR r - 1 rO H 

620 2(0) = /(0) TA< ! ):K/( i ) 

625 NLXr i 

630 FOR r---f1 N! ! iU M 
635 Zi'O) -Z(0) HU I ^ 

o40 NL-xr 1 

650 r ^0)-==F"( I )-TO!TL^ ! ) 

660 AINi“(5y!l) 

670 IF H -2 f)0 rn 700 

680 GO TO 660 

700 Tlhi-(Rn 

BOO i!" W i i GO TO 700 

820 GO TO 500' 

830 STOP 

1200 REM GET Hit-ilRrNCE lOOATlOi] 

1205 C(0> -I 

1210 C(l)-0 

J220 FOR [.“1 TO 0 

1230 Ca)~r<J )~A( L) 

J240 NEXT L 


^OOR QUALirp 
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■*1245 IF M~1=0 00 TO 1600 
1250 C(2)==0 
J260 FOR I=-i TO M~1 
1270 FOR J=I41 TO M 
1280 C(2)-C(2)-}Aa))i<A< J) 
1290 NEXT J 
J300 NEXT I 

1310 -IF H-2“0 GO TO 1600 
1320 C(3)=-0 




]330 FOR I--1 TO N-2 
L340 FOR J=-C-H 10 M-t 
1350 FOR Il-J+L ro M 

J360 C(3)=-C{3)-A(i)^A( J)>KA( I 1) ) 

d362 C(3>-C<3)4Aa)KU( J>KD<1J )+B(] .Oy.'IK I I) 
J370 NEXT II 
1380 NEXT J 
1390 NEXT I 


3 400 IF H-3-=-0 GO TO 1600 
I410 C(4)=0 
1420 FOR [=-| ro H-3 
1430 FDR J=l+i 10 h 2 
1440 FOR I1==J-U TO M-J 
1450 FOR JL~[l+j TO M 

1460 C(4)==C(4)+Aa)>!cA( J)^A< J l)JKA( Jl) -B ( I ) ( J ) iKA ( 1 1 ) JKA ( J J ) 

1461 C(4)==C(4)”A( M)KB( J)-T{B( J l)-^<AkJl)-B(Il5HA< IliAC JL) 

3 462 C(4>=C(4)-A(1>-KA( J3^6< mtBf Jl)lB<.I>)i<B( J3;KB( IJ >>!cDCJ! ; 
1463 C ( 4 ) --C ( 4 ) -A ( [J )KB ( J ’> )KA ( L 3 ) I. J3 J -B ( I J iHA ( J ) )!<A ( ] J ) >!tB ( J I ) 
1470 NEXT J1 

1480 NLXT 11 
1490 NEXr J 
3500 NEXr I 

1510 IF M"4=0 00 ro 1600 


I 520 C ( 5 )-- -ft a )-KA< 2 ) >KA ( 3 ) ( 4 ) >KA ( 5 3 -f B U 3 >!<B (2 ) 5icA C 3 J >KA ( 45 *A (5 


1521 C(5 5=-C(53-fA(l 3)KB(2 3 ■i<B^v3KAc4 3 4^A>v5) 
3 522 C(5)=-C(5)+Ba)'iiA(2) };B(3))}:A( 4 5'1'Ak':.'' 
1523 C(5j-C(G5+A(J 5 IrA ( 2 3 < 3 '> KlU 4 ) >ivA ( 5 '> 

1529 C(5 3-=C(5)-B( L)5KB^2'-)KBk3''iB^ 4)KAk5^ 
1 53b C V 5 ^ =C ( 5 ) f A (1 5 < 2 5 >rA ( 3 5 tX-i ( 4 '• ( 5 3 

3 540 C(5)=^CC5)+BU )TA<2>4<A^3)!!^E^4) KAkG) 
3600 

1610 FDR T=-0 TO M 
1620 M-KL+G<I'‘ 

1640 NEXT 1 
3 650 REHJRN 
3 700 LNli 


READY 


125 



1. Report No. 2. Government Accession No. 

NASA CR-145316 


4. Title and Subtitle 

APPLICATION OF DIGITAL CONTROL TO A MAGNETIC MODEL 
SUSPENSION AND BALANCE MODEL 


3. Recipient's Catalog No. 


5. Report Date 

January 1978 


6 Performing Organization Code 


7. Author(s} 

P. B. Luh, E. E. Covert, H. P. Whitaker, S C. W. Haldeman 


9. Performing Organization Name and Address 

Massachusetts Institute of Technology 
Aerophysics Laboratory 
Cambridge, MA 02139 


12 Sponsoring Agency Name and Address 

National Aeronautics & Space Administration 
Washington, DC 20546 


8 Performing Organization Report No. 

MIT TR 198 


10 Work Unit No. 


11 Contract or Grant No. 

NSG-1292 


13. Type of Report and Period Covered 

Final Contractor Report 
4/1/76 - 11/30/77 


14. Sponsoring Agency Code 



16. Abstract 

The feasibility of using a digital computer for performing the automatic control 
functions for a magnetic suspension and balance system for use with wind tunnel models 
is investigated. Modeling is done using both the NASA-MIT prototype MSBS and a one- 
dimensional magnetic balance. A microcomputer using the Intel 8080 Microprocessor is 
described and results are given using this microprocessor to control the one-dimension- 
al balance. 

Hybrid simulations for one degree of freedom of the MSBS were also performed and 
are reported. It is concluded that use of a digital computer to control the MSBS is 
eminently feasible and should extend both the accuracy and utility of the system. 


17. Key Wordj (Suggested by Author($)] 

Magnetic balance 

Digital control 
Wind tunnel testing 


18 Distribution Statement 


Unclassified - Unlimited 


19. Security Oassif {of this report) 

Unclassified 


20 Security Dassif (of this page) 

Unclassified 


21 No of Pages 

124 


22 Price* 

$6.50 


» 

For sale by the National Technical Information Service. Springfield. Virginia 22151 













